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ABSTRACT 

We measure the sky-projected stellar obliquities (A) in the multiple-transiting planetary systems KOI-94 and 
Kepler-25, using the Rossiter-McLaughlin effect. In both cases the host stars are well-aligned with the orbital 
planes of the planets: for KOI-94 we find A = -14± 13°, confirming a recent result by Hirano and coworkers; 
and for Kepler-25 we find A = 2 ± 5°. 

There are now a total of five obliquity measurements for host stars of systems of multiple transiting plan- 
ets, all of which are consistent with spin-orbit alignment. This alignment is unlikely to be the result of tidal 
interactions, because of the relatively large orbital distances and low planetary masses in the systems. In this 
respect the multiplanet host stars differ from hot- Jupiter host stars, which commonly have large spin-orbit mis- 
alignments whenever tidal interactions are weak. In particular the weak-tide subset of hot-Jupiter hosts have 
obliquities consistent with an isotropic distribution (p = 0.6), but the multiplanet hosts are incompatible with 
such a distribution (p ^ 10"^). This suggests that high obliquities are confined to hot-Jupiter systems, and 
provides further evidence that hot Jupiter formation involves processes that tilt the planetary orbit. 
Subject headings: techniques: spectroscopic — stars: rotation — planetary systems — planets and satellites: 
formation — planet-star interactions — stars: individual (Kepler-25) — stars: individual 
(KOI-94) 



1. INTRODUCTION 

Over the last few years, many stars with exoplanets have 
been found to have high obliquities, i.e., large angles be- 
tween the stellar equator and the planet's orbital plane (e.g. 
Hebrard et al. 2008; Winn et al. 2009; Queloz et al. 2010; 
Collier Cameron et al. 2010b; Moutou et al. 2011; Johnson 
et al. 201 1 ; Albrecht et al. 2012b). However, for practical rea- 
sons almost all of the measurements have been made for stars 
with hot Jupiters. Systems with smaller planets, longer-period 
planets, and multiple planets remain relatively unexplored. 

For the hot Jupiters, Winn et al. (2010) and Albrecht et al. 
(2012c) found evidence that the obliquities of many systems 
have been affected by tidal evolution: the systems for which 
one would expect planet-star tidal interactions to be strongest 
are preferentially found to have low obliquities. Systems 
with weaker tides have a more random obliquity distribution. 
This suggests that at the time of hot Jupiter formation, be- 
fore tides had any opportunity to act, the orbital planes were 
only loosely correlated with the equatorial planes of their host 
stars. This in turn provides evidence that whatever "migra- 
tion" process produces hot Jupiters also causes their orbits to 
be tilted away from the initial plane of formation, favoring 
scenarios such as planet-planet scattering or the Kozai effect 
over the once dominant paradigm of gradual inspiral within 
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the protoplanetary disk. 

The interpretation of the hot-Jupiter results is not settled, 
though, because the possibility remains that high obliquities 
are a generic feature of planetary systems, not specific to hot 
Jupiter migration. There are several proposed mechanisms 
for tilting a star relative to its protoplanetary disk: chaotic star 
formation (e.g. Bate et al. 2010; Thies et al. 2011), magnetic 
star-disk interactions (Lai et al. 2011; Foucart & Lai 2011), 
torques due to internal gravity waves (Rogers et al. 2012), 
and torques due to neighboring stars Batygin (2012). In these 
scenarios, we should observe high obliquities not only in hot- 
Jupiter systems but also in a broader class of planetary sys- 
tems. 

One may test this idea by measuring stellar obliquities in 
systems with multiple transiting planets. In such systems the 
planets' orbits are likely to be coplanar, and presumably mark 
the plane of the protoplanetary disk out of which the planets 
originally formed. If these systems have low stellar obliqui- 
ties as a rule, then the high obliquities in hot- Jupiter systems 
are probably due to planet migration. If instead the obliquity 
distribution of multiple-transiting systems is similar to that of 
hot-Jupiter systems, then the obliquities are clues to more gen- 
eral processes in star and planet formation, and not specific to 
hot Jupiters. 

The first multiple-transiting system for which the projected 
obliquity was measured was Kepler-30 (Sanchis-Ojeda et al. 
2012). The authors took advantage of Kepler photometry 
and the occurrence of star-spots to measure the projected 
obliquity. More recently Hirano et al. (2012b) measured the 
projected obliquity in KOI-94 making use of the Rossiter- 
McLaughlin (RM) effect, and Chaplin et al. (2013) con- 
strained the obliquities of Kepler-50 and Kepler-65 using as- 
teroseismology. AH of these systems were found to be consis- 
tent with good spin-orbit alignment. 

In this work we present an obliquity determination for the 
Kepler-25 multiple-transiting system, as well as an indepen- 
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TABLE 1 

Relative radial-velocity measurements 



TABLE 2 
Characteristics of KOI-94 



System 


Time 


RV 


Unc. 




(BJDtdb) 


(ms-') 


(m s-') 


Kepler-25 


245576 L775 13 


-7.07 


3.53 


Kepler-25 


2455761.78302 


-5.47 


3.45 


Kepler-25 


2455761.85205 


-6.80 


3.77 



KIC 

Kepler magnitude 

Tcff 

Metallically, [Fe/H] 

Projected stellar rotation speed, vsini* 

Macroturbulence parameter, 

Stellar radius, Rp 

Stellar mass, Mn 



6462863 
12.2 

6182(58) K 
4.182(0.066) 
0.0228(20) 
7.3(0.5) km S-' 
4.65 kms"' 
1.52(14) Rq 
1.277(50) Mq 



dent observation of the KOI-94 system. Between these and 
the previously pubHshed measurements, there are now five 
muhiple-exoplanet systems for which we have information 
about the stellar obliquity (and of course the Solar system 
provides a sixth multiple-planet system). We are now in a po- 
sition to make a statistical comparison between the multiple- 
planet systems and the hot- Jupiter systems. 

The plan of this paper is as follows. The observations are 
described in Section 2. Section 3 presents the results for KOI- 
94. Section 4 gives the results for Kepler-25, which were 
obtained with two independent analysis methods because of 
the relatively challenging nature of the detection. The first 
method involved analyzing the "anomalous radial velocity" 
due to the RM effect (Section 4 1), and the second method 
involved direct modeling of the line-profile distortions (Sec- 
tion 4.2.4). As a test of the latter method, an analysis of 
archival spectra of the HAT-P-2 system is also presented. Fi- 
nally, Section 5 presents statistical comparisons of the stel- 
lar obliquities in multiple-planet systems and hot-Jupiter sys- 
tems. 

2. OBSERVATIONS 

All the observations analyzed here were obtained with 
the Keck I telescope and its High Resolution Spectrograph 
(HIRES; Vogt et al. 1994). We observed KOI-94 during 
the night of 2012 August 9/10, when it was transited by 
KOT94.01. We observed Kepler-25 on two nights coincid- 
ing with transits of the largest planet c (2011 July 18/19 and 
2012 May 31 /June 1). We determined relative radial velocities 
in the usual way for HIRES, by analyzing the stellar spectra 
filtered through an iodine cell. The iodine absorption lines 
cover the wavelength range from about 500 to 600 run. The 
analysis was performed with a descendant of the original code 
by Butler et al. (1996). The RVs of the Kepler-25 and KOI-94 
systems are presented in Table 1. 

3. KOI-94 

The KOI-94 system has been studied in detail by Weiss 
et al. (2013). It harbors four planets in orbit around a late 
F-stai- (Table 2). KOI-94.01 is the largest of these planets, 
blocking nearly 0.8% of the starlight during transits. By an- 
alyzing the light-curve of a mutual planet-planet eclipse in 
front of the stellar disk, Hirano et al. (2012b) showed that the 
mutual inclination between KOI-94.01 and KOI-94.03 is low 
(1.15 ± 0.55°). This suggests that the planets have not been 
dynamically disrupted and that their orbits are faithful indica- 
tors of the plane of their formation. 

By observing the RM effect with Subaru and its High Dis- 
persion Spectrograph, Hirano et al. (2012b) measured a low 
projected obliquity (A = -6;}3°) for the host star, relative to 
the orbit of KOI-94.01. Here, we present an analysis of the 



-A- Data from MAST archive http : / / archive . st sci . edu/kepler / 
t Data from Weiss et al. (2013) 

I Obtained using the SME package (Valenti & Piskunov 1996) 

very same transit, based on data obtained with a different tele- 
scope. To analyze the RVs we used the approach of Albrecht 
et al. (2012c), which is based on the analytic description of 
the RM effect by Hirano et al. (2011) to model the RM ef- 
fect. The model takes into account the rotation of the star, 
macroturbulence (e.g. Gray 2005), thermal broadening, and 
line-broadening due to the finite resolution of the spectro- 
graph. We added the prescription for the convective blueshift 
developed by Shporer & Brown (2011), as implemented by 
Albrecht et al. (2012b). We assumed the convective blueshift 
to have an amplitude of 1 kms"', larger than that of the Sun 
(0.5 km s~') because it is thought that slightly hotter stars such 
as KOI-94 have stronger convective blueshifts (see Shporer & 
Brown 2011, and references therein). 

Along with the RVs, we analyzed the transit photometry 
obtained with the Kepler telescope in its short cadence mode 
(one-minute sampling). We fitted the photometric data with 
the Mandel & Agol (2002) transit model to determine the or- 
bital period (P), the time of a particular mid-transit (T^), and 
the geometric parameters of the transit. The geometric pa- 
rameters are the stellar radius in units of the orbital distance 
{Ri,/a), the planet-to-star radius ratio (Rp/R^,), and the cosine 
of the orbital inclination (cos/q). We assumed a quadratic 
limb-darkening law and allowed both of the coefficients to 
be free parameters. The light curve was fitted simultaneously 
with the RVs. To calculate the anomalous RV due to the RM 
effect, we also assumed a quadratic limb-darkening law, with 
priors on the coefficients based on the tabulated values of 
Claret (2000): mi,rm = 0.35 and M2,rm = 0.35. We allowed 
(mi,rm + M2,rm) to be adjusted, with prior constraint centered 
on the tabulated value of 0.7 and with a Gaussian width of 0. 1 . 
The difference between the two parameters (mi.rm - M2,rm) 
was held fixed at the tabulated value of 0.0 since the differ- 
ence is only weakly constrained by the data and in turn has 
little effect on the other parameters. 

Our prior on the projected stellar rotation speed (vsin/*) 
was based on the determination by Weiss et al. (2013) us- 
ing Spectroscopy Made Easy (SME; see Table 2). This spec- 
troscopic modeling code is described by Valenti & Piskunov 
(1996); Valenti & Fischer (2005). Our confidence inter- 
val for this prior was enlarged to 1.5 kms"', because the 
vsin/* value measured from the broadening of stellar absorp- 
tion lines might not be fully representative for the projected 
rotation speed of the stellar surface area covered during tran- 
sit. For example depending on the impact parameter, differ- 
ential rotation might lead to such a mismatch. Our prior on 
the macroturbulent velocity is based on work by Gray (1984). 
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Fig. 1. — Rossiter-McLaughlin effect for KOI-94. Left: Apparent RV variation spanning the transit of planet KOI-94.01 on 2012 August 9/10. Black error 
bars show internal uncertainties as estimated by the RV measurement routine; gray error bars show uncertainties after adding a 'stellar jitter' term in quadrature 
with the internal uncertainties to obtain a reduced of 1- In th£ upper left panel, the solid curve represents the best-fitting model. The lower left panel displays 
the residuals, with the light and dark gray bars indicating the calculated times of first, second, third and fourth contact. Right panels: Parameter distributions, 
based on an MCMC analysis. The main plot shows the posterior in the v sin u -A plane, with dark shading indicating high probability. The black, dark gray, and 
light gray lines outline the two-dimensional 68.3 %, 95 %, and 99.73 % confidence limits. The one-dimensional projections of this posterior are shown on the 
upper and right sides. We find A = -14± 13° and I'sini* = 7.3 it 0.6 km s"'. The dashed line indicates the prior knowledge on I'sinr* (7.3± 1.5kms-') which 
was adopted for this analysis. 



From their equation (4) we obtained a macroturbulence ve- 
locity (0 of 4.9 kms"', for which we assumed a confidence 
interval of 1 .5 km s"' . 

To isolate the RM signal one must subtract (or model si- 
multaneously) the orbital RV variation. One possibility is to 
subtract a model of the orbital RV variation based on the RV 
semiamplitude (/T*) obtained by Weiss et al. (2013), which 
was based on RVs obtained sporadically over several months. 
We did not choose this approach out of concern that apparent 
RV variation due to starspots or other astrophysical noise can 
depend strongly on timescale. Starspot-induced signals, for 
example, can introduce slow drifts in the measured RV signal 
on a particular night. Such a signal would be averaged out 
in a data set obtained over many stellar rotation periods. See 
Albrecht et al. (2012c) for examples of this effect and how it 
can influence measurements of A. Consequently we did not 
apply a prior constraint on K^, or on the systemic velocity 7 in 
our analysis.' 

For the estimation of the uncertainty intervals, we used an 
MCMC algorithm (see, e.g., Tegmark et al. 2004). In Table 
3 we report the results derived from the posterior, where the 
quoted uncertainties exclude 15.85 % of all values at both ex- 
tremes, and encompass 68.3 % of the total probability. Fig- 
ure 1 presents the measured RVs, the best fitting model and 
the posterior in the vsin/,t-A plane. The key result is A = 
-14 ± 13°, a low projected obliquity between the orbital an- 
gular momentum of KOI-94.01 and the angular momentum of 
the stellar rotation. Hirano et al. (2012b) found A = -6;!;j[°, 
which is consistent with our results. However, there were 
some differences in the methods of analysis. Hirano et al. 
(2012b) did not account for correlations between the uncer- 

' See also Isaacson & Fischer (2010) for a discussion of stellar jitter and 
its influence on RV measurements. 



tainty in A and the uncertainties in Tc, Ri,/a, R^/Ri,, /q, or 
they did not model the convective blueshift; and they used RV 
observations obtained on different nights. For a fairer com- 
parison we repeated the analysis of their RVs with the same 
constraints we applied to our own data. In this way we found 
A = -9 ± 14° based on the Subaru dataset. We therefore con- 
clude that two independent data sets support the finding of a 
low stellar obliquity in the KOI-94 system. 

4. KEPLER-25 

The transiting objects in the Kepler-25 system were con- 
firmed to be planets by Steffen et al. (2012), through the detec- 
tion and interpretation of transit timing variations (TTV). The 
system was also recently analyzed by Lithwick et al. (2012), 
who measured masses of 7.13±2.5M0 and 13.1 ±2.6M0 
for the two transiting planets b and c. (For comparison the 
mass of Neptune is 17.15M0.) Table -l gives the basic system 
parameters. Detection of the RM effect for this system is par- 
ticularly challenging because the largest planet c blocks only 
0.13% of the starlight. This leads to a low signal-to-noise 
ratio (S/N) detection of the RM effect. To gain more con- 
fidence in the results we employed two different techniques 
for measuring A, relying on two different portions of the stel- 
lar spectrum. These two measurements are largely indepen- 
dent, although for both methods we use the Kepler photom- 
etry as supporting data. In the first technique (section 4.1), 
we analyze the RV time series which is derived from the io- 
dine region of the spectrum. In the second technique (sec- 
tion 4.2.4) we do not analyze RVs; rather, we directly model 
the deformation of the stellar absorption lines. For this we 
use a method which we originally developed for the analy- 
sis of mutual events in eclipsing star systems (Albrecht et al. 
2007, 2009, 201 la, 2012a). A similar approach has also been 
used by Coflier Cameron et al. (2010a,b); Mifler et al. (2010); 
Gandolfi et al. (2012) and Brown et al. (2012). 
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TABLE 3 

Parameters of the KOI-94 system 



TABLE 4 
KOI-94: Stellar characteristics 



Parameter 


Values 


Parameters mainly derived from photometry 


Mid-transit time Ti [BJDtdb-2400000] 


54965.74127 ± 0.00022 


Period. P [days] 


22.3429513 ± 0.0000092 


Cosine orbital inclination KOl-94.01, cosi'o 


0.0105 ± 0.0019 


Fractional stellar radius, 


0.03786 ± 0.00051 


Fractional planetary radius, Rp/Ri, 


0.06968 ± 0.00035 


Ul +U2 


0.533 ± 0.024 




0.272 ± 0.071 


Parameters mainly deiived from RVs 


Velocity offset, 7 [m s"' ] 


2.5 ± 2.1 


Orbital semi-amplitude, [ms"'] 


80 ± 61 


\/vsin4sinA [v^kms^'] 


-0.64 ± 0.55 


Vvsin/jjCos A [\/kms"'] 


2.59 ± 0.16 


Macro turbulence parameter, ^ [kms"'] 


4.9 ± 1.5 


"l,RM +«2,RM 


0.70 ± 0.09 


Indirectly derived parameters 


Orbital inclination KOI-94.01, lo [°] 


89.40 ±0.11 


Full duration, T14 [hr] 


6.679 ± 0.011 


Ingress or egress duration, T12 [hr] 


0.47 ± 0.015 


Projected stellar rotation speed, vsini* [kms"'] 


7.3 ± 0.7 


Projected spin-orbit angle, A [° ] 


-14 ± 13 



4.L Analyzing the RVs 

The analyis of the Kepler-25 RV time series was similar 
to the analysis of the KOI-94 RV time series. Because for 
Keopler-25 we had RV measurements from two different tran- 
sit nights, separated by nearly one year, we introduced for 
each night a different velocity offset (7) and a different pa- 
rameter to fit the out-of-transit velocity slope (iCt). We used 
a prior on vsini^. based on an SME analysis (see Table 4). 
For C, we used 4.8 ±1.5 kms"', which was obtained in the 
same way as for KOI-94. From the tables of Claret (2000) we 
obtained prior information on the limb-darkening coefficients 
(mlrm = 0.33 and M2,rm = 0.36). As the expected RM signal 
has an amplitude of only a few m s~' , the convective blueshift 
(CB) might have a significant influence on the observed RM 
signal (Shporer & Brown 2011). Thus we allowed the CB 
parameter to vary, with only a weak prior of 1 ±0.5 kms"' 
instead of keeping it fixed as we did for KOI-94. 

Because the planet shows TTVs we used a slightly differ- 
ent approach for incorporating the Kepler photometry into 
our analysis. First we examined the Kepler photometry for 
the two particular transits observed with HIRES, and derived 
midtransit times. From these we computed the ephemeris 
7; = 2455762.03086 ± 0.00050 BID and P = 12!'7203424 ± 
0?00003 days. Second, we measured midtransit times for all 
the Kepler transits and used these to create a single, phase- 
folded, high-S/N transit light curve for Kepler-25 c. This light 
curve was then fitted together with the RVs, providing tight 
constraints on the geometric transit parameters. 

Figure 2 shows the RVs from both nights, and the results for 
vsin/v, and A based on an MCMC analysis of these RVs. The 
results are given in Table 5, second column. Our measurement 
of A = 5 ± 8° is consistent with alignment between the orbital 



KOI 244 

KIC 4349452 

Kepler magnitude 10.73 mag 

Teff 6301(82)K 

logg 4.02(0.1) 

Metallically, [FeH] -0. 10(4) 
Projected stellar rotation speed, vsin/* kms"' 9.5(0.5) 

Stellar radius, Rp 1.36(13)^0 

Stellar mass. My, 1 .22(6)*/© 



*: Data from MAST archive http : //archive . stsci . edu/kepler/ 
t Obtained using the SME package Valenti & Piskunov (1996) 
X Data from Steffen et al. (2012) 

plane of planet c and the stellar equator 

Because the amplitude of the RM effect is « 4 m s~', com- 
parable to typical uncertainty of a single RV measurement, 
we must be skeptical and examine this result further. Why 
does our algorithm find such a small uncertainty interval for 
A, given that the RM signal is so difficult to discern in the time 
series of the RV measurements (Figure 2)? 

We have a great deal of prior knowledge of all the system 
parameters relevant for the RM effect, except for A, allow- 
ing us to predict accurately the expected characteristics of the 
RM signal as a function of A. To first order the amplitude 
of the RM effect is proportional to the covered surface area 
and the projected rotation speed. (See Gaudi & Winn 2007; 
Albrecht et al. 2011b, for a more detailed discussion.) How- 
ever the amplitude of the RM effect also depends on A itself. 
The amplitude of the RM signal is nearly twice as large for 
A = ±90° as it is for A near 0° or 180° (Figure 3). Because of 
this and because there is a hint of a prograde signal in the RVs 
(Figure 2) the low projected obliquity is favored. 

In more detail: we know from the Kepler photometry that 
planet c has a high impact parameter, i.e., it travels near the 
stellar limb throughout the transit. We also know a priori 
that the star has a substantial projected rotation speed from 
the SME analysis. Combining these two pieces of informa- 
tion we know there is no way to make the RM effect vanish.- 
Figure 2 illustrates that we did not make a high-S/N detection 
of the RM effect. We might therefore ask which values of A 
would lead to the smallest RM amplitude, or an RM signal 
which would be easiest to hide by adjusting other parameters 
in our model. We have just noted that the maximum ampli- 
tude of the RM effect is larger for |A| near 90° than for A near 
0° or 180° (Figure 3). This is because limb darkening attenu- 
ates the signal by a factor of 0.6-0.7 for a star like Kepler-25. 
The RM signal has maxima when the planet crosses the stel- 
lar limb where limb darkening is strongest, and the signal is 
zero near mid-transit where limb darkening is weakest. For 
|A| near 90 degrees the maximum RM signal occurs nearer to 
the center of the stellar disk. In addition, the parameters 7 and 
Ki, (the out-of-transit slope) are most strongly covariant with 
A when A is near 0° or 180° (Albrecht et al. 2011b). This is 
because of the time-antisymmetry of the RM signal in such 
cases. In contrast, for A near ±90° the RM signal is time- 

^ If the projected stellar rotation speed were not known, or if the planet had 
a low impact parameter, then it would be possible to reduce the amplitude of 
the RM signal to arbitrarily small values Albrecht et al. (201 lb). 
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Fig. 3. — RM signal for different obliquities. The plot shows the ex- 
pected RM signal for a system like Kepler-25 during transit of planet c. The 
solid line shows the signal for the parameters printed in the second column of 
Table 5 but now with A = 0° . The dashed and dashed dotted lines show the 
signals for A = 90° and A = —90° . The maximum amplitude of the aligned 
signal is smaller than for the strongly misaligned cases. 

symmetric. Finally, differential rotation would also weaken 
the RM signal for A = 0° or 180° compared to the A = ±90° 
signal, though this is a comparatively minor effect 10%). 
Together these factors make it easier to hide an RM signal 
with A = (or 180) than A = ±90°. Therefore, it is possible to 
infer that A must be near 0° or 180° with a sufficiently con- 
straining upper limit on the amplitude of the RM effect. Here 
the data prefer 0° over 180° as there is a hint of a prograde 
RM signal in the data. As mentioned above there is only one 
parameter that is not at least partly constrained by photometry 
or prior knowledge, which is the projected obliquity. This is 
the qualitative explanation for the MCMC result of A = 5 ± 8 
degrees. Note that the a priori knowledge on vsin/v, is crucial 
in this situation (e.g Albrecht et al. 201 lb). 

However it is not completely satisfactory to argue for a low 
obliquity based on the absence of a clear signal. Therefore we 



sought an independent method to detect the RM effect. The 
next section describes our analysis of a wavelength region in 
which no iodine lines are present, and which therefore was not 
used in the determination of the RVs. We used this spectral 
range to make a second independent measurement of A. 

4.2. Measurement of the planet's Doppler shadow 

In this method of analysis we did not use RVs as prox- 
ies for the distortions in the stellar absorption lines. Instead 
we analyzed the line shapes directly to infer A. The transit- 
ing planet selectively distorts the line profile, blocking a cer- 
tain range of velocity components that ordinarily contribute 
to the overall line broadening, a phenomenon referred to as 
the "Doppler shadow" by Collier Cameron et al. (2010a). To 
model this phenomenon we used a code developed for double- 
lined eclipsing binaries (e.g. Albrecht et al. 2007, 2012a). For 
binary star systems, the difficulty lays in the additional set 
of stellar absorption lines originating from the eclipsing fore- 
ground star. In the case of planetary transits, this particular 
difficulty is eliminated, as the planet's emitted light is neg- 
ligible; rather, the challenge stems from the small amplitude 
of the RM signal. In the case of a Jupiter-sized planet such 
as HAT-P-2b, about 1 % of the light is blocked from view. In 
the case of the transit of planet c in front of Kepler-25 only 
0.13 % of the light is blocked. For this reason we will present 
the application of our code to the Jupiter-sized planet in the 
bright (my = 8.7) HAT-P-2 system, as a test case. We then 
proceed to the more challenging Kepler-25 system. 

4.2.1. The method 

We used a two-step process to measure A from stellar spec- 
tra. In the first step we combined the signals from many stellar 
absorption lines into one high-S/N absorption line, which we 
will call the "kernel." The spectrum is modeled as the convo- 
lution of the kernel and a series of (5-functions at the central 
wavelengths of the absorption lines. This was done for each 
individual spectrum. In the second step, we analyzed the dis- 
tortions that are seen in the kernels, which are caused by the 
transit of the planet over the rotating photosphere. As the first 
step is different from the approach used by Albrecht et al. 
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TABLE 5 

Results for the Kepler-25 system. 



Parameter 


Values RV 


Values Distortion 


Parameters mainly derived from photometry 


Mid-transit time Tc [BJDtdb-2400000] 
Period, P [days] 

Cosine orbital inclination KOI-94.01, cosi'o 
Fractional stellar radius, R^/a 
Fractional planetary radius, Rp/Rt 

Ml + U2 
Ul —U2 


55762.0309 ± 0.0005" 
12.72034 ± 0.00003" 
0.0475 ± 0.0008 
0.0540 ± 0.0007 
0.0360 ± 0.0005 
0.560 ± 0.020 
-0.14 ± 0.5 


55762.0308 ± 0.0005" 
12.7203424 (fixed) 
0.0476 ± 0.0008 
0.0540 ± 0.0007 
0.0362 ± 0.0007 
0.560 ± 0.019 
0.04 ± 0.6 


Parameters mainly derived from spectroscopy 


Velocity offset 201 1, 7 [m s"' ] 

Velocity offset 2012, 7 [m s"' ] 

Orbital semi-amplitude 2011, K* [ms"'] 

Orbital semi-amplitude 2012, K* [ms"'] 

V I' sin U sin A [Vkm s^' ] 

V V sin /j, cos A [Vkm s"' ] 

Macro turbulence parameter, ( [kms"'] 

Convective blueshift [kms"'] 

Point Spread Function, PSF [km s"' ] 


-3.7 ± 1.3 

2.7 ± 0.7 
-10 ± 23 
-32 ± 24 

-0.35 ± 0.39 
2.9 ± 0.23 

4.8 ± 1.5 
-1.0 ± 0.5 
0.70 ± 0.10 

3 (fixed) 


b 
b 
b 
b 
b 
b 

4.1 ± 1.1 
-0.8 ± 0.07 
0.65 ± 0.09 

4.2 ± 0.4 


Indirectly derived parameters 


Orbital inclination KOI-94.01, i'o [°] 

Full duration, T[4 [hr] 

Ingress or egress duration, T12 [hr] 

Projected stellar rotation speed, v sin r'^. [km s"' ] 

Projected spin-orbit angle, A [° ] 


87.28 ± 0.05 
2.860 ± 0.009 
0.405 ± 0.014 
8.7 ± 1.3 
7 ± 8 


87.27 ± 0.05 
2.861 ± 0.009 
0.410 ± 0.015 

8.2 ± 0.2'^ 
-0.5 ± 5.7'^ 



Notes — 

" We used the priors P = 1 2:^7203424 ± 0:^00003 and 7i = 2455762.03086 ± 0.00050 BJD, 
as determined using the Kepler photometry of the appropriate transits (see section 4. 1 ). 
Was not determined. 

Here we step directly in v sin u and A as they are less correlated than for RV measurements. 



(2007), we discuss it in detail in the following subsection. 

4.2.2. Measuring high S/N ratio stellar rotation kernels 

Preparing the spectra. — The new algorithm works as follows. 
First we normalize the spectra using observations of fast ro- 
tating B-stars. Specifically we use polynomials fitted to the 
same and adjacent orders in the B-star spectrum to normalize 
the spectral orders of our science target. We use adjacent or- 
ders to remove the influence of shallow spectral lines present 
in the B star spectrum. All spectra are shifted according to the 
barycentric correction, and a correction term derived from the 
measurements of deep telluric lines in the red wavelength arm 
of HIRES. These corrections line up the spectra with an accu- 
racy of 100-200 m s"' (but see also below). Next, a high-S/N 
spectrum is created by averaging over all out-of-transit ob- 
servations obtained during the night. Now each spectrum is 
compared to this high-S/N spectrum to mark and discard bad 
pixels. We also make one final small differential correction 

^ Previously, Albrecht et al. (2007) used the broadening function devel- 
oped by Rucinski (1999). For comparison, we have reanalyzed those data 
using the algorithm presented here, finding equivalent results for A and the 
other system parameters. However, for Kepler-25 we found it easier to create 
a template spectrum using the new approach presented here, because of the 
availability of higher-S/N spectra. The HIRES spectra of Kepler-25 have a 
S/N between 50 and 100 in the wavelength range 398 to 479 nm. 



in the normalization of all spectra. For this we compare all 
spectra to the high-S/N spectrum and fit a third-order polyno- 
mial to the residuals, in which no absorption lines are present. 
Such a polynomial is created for each order and spectrum and 
is used for the correction. 

Creating and refining a template. — To combine all the infor- 
mation contained in the different absorption lines, we need a 
sharp-lined template spectrum matching that of the target star. 
Our approach to obtain such a template — which matches the 
observed spectrum after convolution with a "master" absorp- 
tion line profile — was inspired by Reiners & Schmitt (2003), 
but see also Donati et al. (1997) and Rucinski (1999) for sim- 
ilar approaches to obtain high-S/N kernels. 

It is not only important that all the lines present in the ob- 
served spectrum are also present in the template; it is also im- 
portant that the line depths in the template match the depths 
in the observed spectrum. We establish the appropriate line 
depths in the template spectrum in the following manner. As 
a starting point for this iterative process we query the Vienna 
Atomic Line Database (VALD; Kupka et al. 1999) for Une 
positions and line depths appropriate for a star with the given 
effective temperature and surface gravity (our inputs in this 
case are based on the SME analysis; see Table 4). We now 
adjust the line depths so that a kernel convolved with the line 
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Fig. 4. — Small portion of the normalized spectrum of Kepler-25. The thin gray line represents a small part of the third observed spectrum during the night 
of 2012 May 31/June 1. The dark line represents the convolution of the template (marked by the dashed lines) with the optimized kernel. The optimized kernel is 
shown in the lower right comer The lower panel shows the difference between the observed spectrum and template convolved with the kernel. For each spectrum 
of Kepler-25, a total of about 900 A were analyzed in the same fashion as shown in this 4 A section. 



list gives the best fit (smallest x^) to the high-S/N spectrum 
we had obtained above. This kernel is created simultaneously 
with the optimization of the line depths. At this stage the 
kernel is purely phenomenological, and is not subject to any 
physical boundary conditions. It consists of a number of free 
parameters. Each parameter represents one pixel along the 
dispersion direction of the spectrograph, translated into ve- 
locity for the wavelength region of interest (« 1.3 kms"' at 
390-480 nm). The number of pixels is chosen such that the 
range of velocities that is covered is about twice as large as the 
vsin/j, of the star For HAT-P-2 we have 61 free parameters 
and for Kepler-25 we have 33 free parameters. To increase 
the speed of the computation each order is split into several 
sections, and for each section the best fitting line depths and 
kernels are found separately. 

After this initial round we create an average kernel out of all 
kernels from all sections in all orders. Here the kernels from 
each section are weighted by the blaze function (for which 
we use the B-star spectra as proxy) of that section, and the 
equivalent width of the absorption lines in that section. We 
exclude regions at the short and long wavelength end of each 
order, regions dominated by deep, large lines, and regions for 
which clearly no good fit was achieved, i.e. lines are missing 
in the template. About 10% of the available spectral range 
was excluded for one or another of these reasons. 

This newly obtained average kernel is now used for all sec- 
tions in a second round, during which only line depths are 
adjusted. In the following round only the section kernels are 
optimized. These last two steps are repeated two more times 
to optimize the line depths. Finally the template with the op- 
timized line depths is saved for later use, while the high-S/N 
spectrum as well as the average kernel are discarded. 

Measuring high S/N ratio kernels. — Now we use the depth- 
optimized template to obtain an average absorption line ker- 
nel in each section of all observations. Figure 4 displays one 
section of one Kepler-25 observation during the transit night. 
Also shown is the kernel derived from this section. Using the 
same weighting scheme as used above we can now create an 
average kernel for each observation. 



4.2.3. Analyzing high-S/N kernels 

The kernels obtained in the last section are then analyzed 
with the same code which was used by Albrecht et al. (2007). 
In short we pixelate the visible stellar surface and assign to 
each pixel a radial velocity, based on contributions from stel- 
lar rotation, macroturbulence, and the convective blueshift. At 
each phase of the transit we integrate over the exposed por- 
tion of the stellar surface to obtain the stellar absorption-line 
kernel. In this step we assume a quadratic limb-darkening 
law. Finally the absorption line is convolved with a Gaus- 
sian function representing both micro-turbulence in the pho- 
tosphere and the point-spread function (PSF) of HIRES. 

Changes in the spectrograph PSF during the transit nights. — Be- 
fore the model absorption lines can be compared to the mea- 
sured kernels one last step has to be taken which is specific 
to these observations. The spectra are obtained with a slit 
spectrograph and we therefore have to take into account pos- 
sible changes in the PSF throughout the night. This is because 
small changes in telescope guiding cause variability in the il- 
lumination of the slit, and consequent changes in the PSF of 
the spectrograph." Changes in illumination of the slit might 
shift the PSF, sharpen or widen it, as well as introduce skew- 
ness and higher frequency terms. As a measured spectrum is 
convolved with the PSF such changes do effect the measured 
absorption lines directly. 

We are interested in a time-varying signal: the distortion 
due to the transiting planet. Therefore it is important to com- 
pensate for changes in the absorption lines due to changes in 
the PSF. It is not crucial to know the PSF itself. In the next 
paragraph we explain how we deal with these potential shifts 
and stretches, and how we deal with higher-order changes by 
interpolating between observations obtained just before and 
after transit. 

To compensate for changes in the absorption lines caused 

^ In principle some information about the time-variable PSF is contained 
in the solutions provided by the RV-measuring code, which is based on a fit 
to the spectra in the iodine wavelength range from 500-600 nm. However we 
did not thoroughly investigate such an approach as the PSF is also expected 
to vary with position on the spectrograph CCD. 
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Fig. 5. — Doppler shadow of HAT-P-2b. During the planetary transit, part of the stellar photosphere is hidden from view, distorting the rotationally-broadened 
stellar absorption lines. The left upper panel (A) shows the time-variable planet shadow during a transit in the HAT-P-2 system. The dashed lines indicate times of 
first and last contact. At the beginning of the transit, parts of the approaching stellar surface area are hidden from view, and therefore blueshifted light is hidden. 
At the end of the transit, redshifted light is hidden. In addition to the transit signal, some artifacts are visible. They originate mainly form hour-to-hour ("fast") 
changes in the wavelength position of the kernels. A mismatch in wavelength between the out-of-transit kernel and the current kernel leads to a deficit on one 
side of the kernel (dark) and a positive residual on the other side (light color). Panel (B) shows the results after correcting for these fast changes, as described 
in the text. Here, for illustration purposes, we only used an average kernel based on the post-transit data, rather than he interpolation between pre-transit and 
post-transit kernels, as described in the text. There is a continuous buildup of difference between the observed and assumed kernels towards the beginning of the 
night ("slow" changes). Panel (C) shows the results when the assumed kernel is based on inteipolation in time between the pre-transit and post-transit kernels. 
Panel (D) shows the residuals after our best fitting model for the planetary transit is subtracted from the data shown in panel (C). There are some low-level patterns 
in the residuals, likely originating form changes in the PSF which are not modeled by our algorithm. To facilitate comparison between the different panels, the 
kernels form each system have been normalized to a height of unity (see Figure (> and Figure 7). The grayscale bar next to each panel indicates the signal strength 
on the same scale. It is interesting to note that the depth of the HAT-P-2b velocity signal is more than twice as deep as the photometric transit signal. For rapidly 
rotating stars the depth of the deformation is not proportional to the square of the radii ratio, but is better approximated by the ratio itself. 



(A) planetary shadow (B) residuals 




velocity (l<m s ') velocity (l<m s ') 

Fig. 6. — Doppler shadow of Kepler-25c. The same as panels (c) and (d) in Figure .'i, but this time for the transit on 2012 May 31/June 1 of Kepler-25c. In the 
left panel (A) one can see the signature of a distortion traveling from negative RV towards positive RV throughout the transit: the signature of a prograde orbit. 
Modeling the distortion gives A = —0.5 it 5.7°, indicating good alignment between the projections of the stellar and orbital spins. Panel (B) shows the residuals 
after subtraction of our best-fitting model. 



Low obliquities in compact multiplanet systems 



9 



average line 




T 1 1 -r 



-20 -10 10 20 

velocity (km s"') 

Fig. 7. — Comparison of overall line shape to model. The solid line 
shows the average line shape of Kepler-25 as measured with the first three 
observation during the transit May/June 31/1, 2012. The gray dashed line 
shows our best-fitting model of this line. The model includes quadratic 
limb-darkening, stellar rotation, micro and macroturbulence, the convective 
blueshift, and a PSF which is modeled as a simple Gaussian. Apart from the 
(jpsF parameter, all other parameters are also simultaneously used to find the 
best-fitting model for the planet shadow (Figure Our results for HAT-P-2 
are not shown here but are of similar quality. 

by changes in the PSF we performed the following steps. We 
take the mean of the first few spectra, and the mean of the last 
few spectra, during a transit night. (If the PSF of the spectro- 
graph would have been stable and our correction for any RV 
changes would have been perfect, then these should be identi- 
cal, assuming no stellar activity). Next we linearly interpolate 
in time between these to create an absorption line appropriate 
for the time of each individual observation. These lines do not 
contain a transit signal, but include slow monotonic changes 
in the PSF. To isolate the transit signal we subtract the mea- 
sured kernels from these interpolated lines. 

To correct for low-order fast changes in the PSF we allow 
the kernel of each observation to shift in velocity space and 
we further allow for a scaling in the velocity scale of the mea- 
sured absorption line. This leads to two free parameters for 
each observation, which are evaluated each time the observa- 
tions are compared to a specific model. The drifts are less then 
200 m s"' and the scaling in velocity space is always less than 
0.5 %. See Figure 3 for an illustration of how these different 
corrections influence the signal. 

With this scheme we can compensate for a constant un- 
known PSF, slow changes in the PSF, as well as fast low-order 
changes in the PSF. We do not attempt to correct for fast high- 
order changes in the PSF as these would likely be correlated 
with the planet transit signal, which is itself a higher-order 
change in the stellar absorption lines. 

In addition to analyzing the changes in the absorption lines, 
we also use the observations taken before transit and compare 
them to the model line. This gives additional constraints on 
vsin/j,, and limb darkening. Here the PSF is modeled as 
a simple Gaussian with only one free parameter, the width 
(cpsf)- No accurate modeling of the PSF is required as we do 
not attempt to identify a small transit signal, but rather sim- 
ply to measure the line width. See Figure 7 for a comparison 
between a measured absorption line and our model. The eval- 
uation of the planet shadow and the out-of-transit line is done 
simultaneously. 



HAT-P-2. — As our technique had to this point only been 
used for binary star systems, we applied it to the HAT-P- 
2 system before using it on the more challenging Kepler-25 
dataset. The spectra of the HAT-P-2 system were obtained 
and analyzed by Winn et al. (2007) and the RVs were re- 
analyzed by Albrecht et al. (2012c). For the line-distortion 
analysis we used the spectral region from 443 to 455 nm. We 
used the same photometric priors as used by Albrecht et al. 
(2012c), where a low obliquity was found (A = 9 ± 10°). With 
our analysis of the spectra we find A < 10° (Figure 5). We 
thereby gained confidence that our algorithm to extract pro- 
jected obliquities directly from modeling stellar absorption 
lines also works for the case of planetary transits observed 
with a slit spectrograph. 

4.2.4. Kepler-25: changes in stellar absorption lines 

Using the scheme described above we then analyzed the 
Kepler-25 spectra obtained during the two transit nights. We 
used the spectral region from 398 to 479 nm, blueward of the 
iodine lines. We did this separately for each transit night. We 
first focused on the transit night 2012 May 3 1/June 1 . We used 
the same composite photometric light curve as used in our 
analysis of the RVs, and left P fixed at the value determined 
from the photometry. With this method of analysis there is 
no need to determine K^, and 7 for the transit nights, as had 
been necessary for our analysis in Section 4. 1 . This reduced 
the number of free parameters by two; on the other hand there 
was a new parameter ctpsf as described above. We did not 
impose any prior constraints on this parameter. To estimate 
the parameters and their uncertainty intervals we again used 
an MCMC algorithm. 

Figures 6 and 7 show the comparison between the data 
and the best-fitting models. Our results for vsin/* and A are 
shown in the left panel of Figure S and printed in column 3 
of Table 5. These are consistent with the results found using 
the RVs of both data sets (Section 4. 1 , Table 3 column 2). 
In particular the results for A are consistent with each other. 
Analyzing RVs obtained with the iodine technique we ob- 
tain A = 5 ± 8°. Analyzing the change in the absorption lines 
in the blue part of the spectrum we obtain A = -0.5 ±5.7°. 
These results are independent from each other as different 
wavelength regions of the obtained spectra have been used 
(although the supporting photometric data was the same in 
both cases). In addition, the two independent measurements 
of vsin/^ are consistent with each other. To illustrate the con- 
sistency of the results we show in the right panel of Figure 8 
the expected RV signal from our best solution to the distortion 
of the absorption lines in the blue spectra. We further show 
the RVs measured during the two transits in the red part of the 
spectra, and which are displayed in Figure 2. We also show 
the difference between the measured RVs and the implied RV 
changes from the distortion. These signals are simply plotted 
on the same axes; they were not adjusted to match each other. 
The agreement of these two different approaches lends addi- 
tional confidence to the conclusion that the projected stellar 
obliquity is low in Kepler-25. 

However, for the first transit night (2011 July 18/19) there 
was no secure detection of the Doppler shadow. What might 
have prevented a detection in this case? The main differ- 
ence between the two datasets is a difference in the range of 
airmasses through which the observations were made (Fig- 
ure 9, upper panel). During the 201 1 transit the airmass in- 
creased from 1.1 to 1.9. In contrast, in 2012 the transit was 
observed at low zenith angles, always below an airmass of 
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Fig. 8. — Kepler-2S results form the analysis of the absorption lines and comparison to the RVs. The left panels show our results from the analysis of the 
distortion in the stellar absorption lines. We measure A = -0.5 it 7° and vsiniV = 8.2 it 0.2 km s"'. To allow a comparison to the RV data, the solid line in the 
right panel shows the anomalous RV signal that is implied by the best-fitting model of the line distortions. It is seen here to be compatible with the RV data even 
though the RV data was not used directly to constrain this model. (For this visual comparison, the out-of-transit RV trends were subtracted from the measured 
RVs.) The lower panel shows the difference between the measurements and the model, illustrating the good agreement between the line-distortion method and 
the anomalous-RV method for characterizing the planetary transit. 
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Fig. 9. — Changes in airmass during the Kepler-25 observations. The 

upper panel shows the changes in air-mass for the observations in 201 1 (filled 
symbols) and 2012 (open symbols). The transit interval is indicated by the 
gray bars. Lower panel: The measured RVs for telluric lines on the red CCD 
of HIRES. 



1.1. The large variation in airmass during the 2011 obser- 
vations strongly increases the difficulty of modeling the ob- 
served spectra. This is because a change in the angle under 
which an object is observed can lead to changes in the slit il- 
lumination which in turn changes the effective PSF. This can 
most readily be seen in the lower panel of Figure 9 where 



the apparent shift in radial velocity of telluric lines is plot- 
ted against time. The apparent changes are about « 1 kms"' 
for the night where we detect no signal and only « 150 m s"' 
where the changes in airmass are low and where the transit 
signal was detected. Telluric lines are intrinsically stable to at 
least a few tens of ms"' (e.g. Gray & Brown 2006; Figueira 
et al. 2010). 

The apparent changes in radial velocity of the telluric lines 
are measurement artifacts, caused by changes in the PSF, 
which also apply to stellar lines. We suspect that on top of 
the RV shifts higher order changes occur, which our algo- 
rithm cannot correct with sufficiently accuracy to allow for 
the detection of the small distortion induced by the planetary 
transit. During the observations of our test system HAT-P-2 
the air-mass also changed significantly, from 1 .2 to 2, and the 
observed RVs of telluric lines changed by w 1 km s"' . The 
residuals of 0.8 % in our test system (Figure 5, panel D) are 
larger than the expected distortion of 0.2 % in the absorption 
lines in Kepler-25 (Figure 6, panel A). The signal of HAT- 
P-2b was nevertheless detected because of its relatively large 
amplitude; but the nondetection of a transit signal during the 
201 1 observations of Kepler-25 is not surprising in this con- 
text. 

5. COMPARISON WITH HOT-JUPITER SYSTEMS 

In this section we will analyze our results for KOl-94 and 
Kepler-25 together with the results for three additional multi- 
transiting systems to try and clarify the interpretation of the 
high obliquities seen in hot- Jupiter systems. The other three 
systems are Kepler-30, 50, and 65. Using the occurrence 
of star spots, Sanchis-Ojeda et al. (2012) measured A to be 
1 ± 10° in Kepler-30. Chaplin et al. (2013) measured the incli- 
nation of the stellar rotation axis towards the observer (/*) for 
two additional systems, via asteroseismology. For Kepler-50 
they found = 82^^° and for Kepler-65 they found = 8 1!?^" . 
An estimate of the stellar inclination in Kepler-50 had also 
been obtained earlier by Hirano et al. (2012a), using a combi- 
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Fig. 10. — Projected obliquity (either A or as a function of the relative tidal-aUgnment timescale, for hot- Jupiter and multiple-planet systems. The 

systems are plotted as a function of a simple metric for the expected timescale for tidal dissipation within the star (see Albrecht et al. 2012c eq. 2-3). Stars which 
have temperatures higher then 6250 K are shown with red filled symbols. Blue open symbols show stars with temperatures lower then 6250 K. Stars which 
measured effective temperature include 6250 K in their l-cr interval are shown by split symbols. Systems for which A was measured are indicated by a circle and 
refer to the left-hand axis. Measurements of u are indicated by a square and refer to the right-hand axis. Systems which harbor multiple planets are given dark 
black borders. The systems with short tidal timescales are seen to be well-aligned. All of the multiple-planet systems are well-aligned despite having weak tidal 
dissipation. 



nation of estimates for the rotation period, stellar radius, and 
vsin/v,. Their result was less constraining than, but compatible 
with, the result of Chaplin et al. (2013). 

High obliquities: planet migration or star-disk evolution? — For 
HJ systems, evidence has accumulated that the stellar obliq- 
uities varied over a very large range when the gas giants ar- 
rived near their host stars (Winn et al. 2010; Albrecht et al. 
2012c). This has been taken as evidence that the orbital plane 
of the planet has changed, presumably via the same mech- 
anism which also changed its orbital distance. However, as 
mentioned in the introduction, there are other mechanisms 
which might create a misalignment between the stellar ori- 
entation and the planetary orbit. In multiple-transiting planet 
systems there is reason to think that the orbits still trace the 
plane of the disk out of which they have formed. There- 
fore measuring the obliquities in these systems lets us learn 
about the degree of alignment between protoplanetary disks 
and stellar spin axes. 

If we find that the distribution of obliquities for multiplanet 
systems is closer to alignment, than the distribution of obliq- 
uities for HJ systems, then this would indicate that the large 
obliquities in HJ systems are caused by the evolution of the 
planets' orbits. If on the other hand we find that the distri- 
bution in obliquities for multiplanet systems is similar to the 
distribution of obliquities for HJs then we the measured obliq- 
uities are not necessarily related to planet migration. 

We note that the host stars in both groups, close in gas gi- 
ants and multiple planet systems, are on the main-sequence 
and cover the spectral classes from F to K. The only read- 



ily apparent difference between these systems is that for one 
group, several planets are found on compact coplanar orbits, 
within in the other group there are solitary transiting gas gi- 
ants. 

The influence of tides. — Before we can compare the two dis- 
tributions of obliquities we first need to know if tides have 
dampened the obliquities. It would be advantageous to only 
include systems which have not undergone any significant 
tidal influence, instead of attempting to model the influence 
of tides on the obliquity. 

To check which systems might be influenced by tides, we 
calculate the same two alignment timescales presented in Al- 
brecht et al. (2012c) for the multiple planet systems. Either 
approach to calculating the timescale leads to the conclusion 
that tides have probably had little or no influence on the stellar 
obliquities in the five multiplanet systems. In Figure 10 we 
show the results for the timescale which is calibrated based 
on binary-star data (Albrecht et al. 2012c). In the same Fig- 
ure we also show all the hot- Jupiter systems from Albrecht 
et al. (2012c), after including some newly published measure- 
ments. We included the new measurements for WASP-32 and 
WASP-38 from Brown et al. (2012), HAT-P-17 Fulton et al. 
(2013), updated the value for CoRoT-11 from Gandolfi et al. 
(2012), and updated the value for WASP- 19 from Tregloan- 
Reed et al. (2013). ^ As mentioned we also added the A and i^, 

^ Recently Hebrard et al. (2013) measured A in the WASP-52 system. We 
do not include this result in the current study. The reasons are similar to the 
reasons for which we excluded some systems in the Albrecht et al. (2012a) 
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measurements for the multiple transiting planet systems. The 
left-hand y-axis indicates A and the right-hand y-axis indicates 
Indication of good alignment is a low value in A or a large 
value of The long tidal realignment timescales in the stud- 
ied multiple planet systems is due to the long orbital periods 
and the small masses of the planets. 

According to Figure 10 it does seem unlikely that any of the 
multiple planet system was influenced by tides. But which hot 
Jupiter systems should be included in the comparison? As we 
do not have a clear cut criterion we will use three samples: (1) 
All hot Jupiters; (2) all hot Jupiters with t equal or larger to 
the T of the first clearly misaligned system (r > 10^'^); and (3) 
only hot Jupiters which have timescales r equal to or larger 
than the shortest t found amongst the multiple-planet systems 
(t= lO^**). 

Comparing the distributions. — In the regime of weak tides, the 
hot Jupiter results appear to be nearly random. Therefore we 
first ask: could either population be drawn from an isotropic 
distribution on a sphere? For hot Jupiters we have only mea- 
surements of A, the projected obliquity. We can therefore 
compare these measurements to a distribution in A for the 
isotropic case using a Kolmogorov-Smirnov (K-S) test (e.g 
Press et al. 1992). For case (1) the K-S test suggests that there 
is negligible probability that the A measurements of all hot- 
Jupiter systems are drawn from an isotropic distribution. For 
case (2) there is still only a 0.04 % probability that the results 
are drawn from an isotropic distribution. However for case (3) 
we find there is a 61 % chance that this distribution is consis- 
tent with an isotropic distribution in A. Figure 1 1 shows the 
cumulative distribution in A for these HJs and the expected 
distribution in A for an isotropic distribution. 

For the multiple planet systems we have two distinct mea- 
sures of obliquity, A and which cannot be translated into 
each other (at least not without already assuming a distribu- 
tion, see Fabrycky & Winn 2009). Therefore we use a Monte 
Carlo approach instead of a K-S test. We create a distribution 
of obliquities which has a uniform distribution in A (Kepler- 
30, KOl-94, Kepler-25) and in cos/,, (Kepler-50, Kepler-65) 
to simulate a isotropic distribution in the obliquities. Form 
these we draw five "measurements" which we compare to the 
three measurements of A and two measurements of i^,. The un- 
certainties in the actual measurements are included as Gaus- 
sian random numbers, every time a comparison is made. In 
particular for the comparison in we use half-Gaussians with 
peaks at 90° and standard deviations derived from the incli- 
nation measurement plus the measurement uncertainty. We 
repeat this experiment 5 x 10^ times. Only in 0.0003% of 
these experiments do we draw sufficiently low values of A, 
and sufficiently high values of to be compatible with the 
measured As and It seems unlikely that obliquities in mul- 
tiple planet systems are drawn from an isotropic distribution. 
A narrow distribution centered near zero obliquity is more ap- 
propriate. We will defer an analysis of which is the exact 
distribution until more obliquity measurements in multiplanet 
systems are available. See Fabrycky & Winn (2009) for pos- 
sible approaches on how a comparison can be made. Such 
an analysis would be interesting as it might shed light on the 
origin of the small (6°) obliquity of the Sun. 

Now combining that 1) multiplanet systems have a differ- 
ent obliquity distribution than systems with single, close-in 
gas giants, 2) planets in multiple planet systems presumably 

study. 
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Fig. 11. — Kolmogorov-Smirnov test on the isotropic distribution of A 
for single, close in, gas giant planets. The solid line shows the cumulative 
fraction of |A| for an isotropic distribution of obliquities on a sphere, all A 
are equally likely a priori. The dashed line shows the cumulative distribution 
for measurements of A in HJ systems. We only included systems with t > 
10^^ (See Figure 10), to avoid systems which are strongly influenced by 
tides. According to the Kolmogorov-Smirnov test there is a 61% chance 
that the projected obliquities of these systems are drawn from an isotropic 
distribution. 

trace with their orbits the plane of the circumstellar disk out 
of which they formed, and 3) we are not able to detect any 
other significant difference between stars which have close 
in giant planets and stars which hosts multiple planets, we 
conclude that the misalignments between stellar rotation and 
planetary orbits are due to changes of the inclinations of the 
orbital planes. Our results disfavor theories which aim to ex- 
plain large obliquities due to a change in the angle between 
protoplanetary disk and the star (e.g. Bate et al. 2010; Thies 
et al. 201 1; Batygin 2012) or changes in the internal structure 
of the star (Rogers et al. 2012). 

Of course it must be acknowledged that a sample of 5 sys- 
tems is not sufficient for a firm conclusion. The systems stud- 
ied here only cover a small parameter space, for example a 
limited range in stellar mass. In other systems, mechanisms 
for tilting stars may be more important. One clue that this is 
indeed the case is the finding that both stars in the Dl Herculis 
system are strongly inclined with respect to their mutual orbit 
(Albrecht et al. 2009). Here however, we have found that the 
evidence to date supports the conclusion that the high obliqui- 
ties of hot-Jupiter systems are due to evolution of the planets' 
orbits. 

Relation to planet migration theories. — If we assume that the 
smaller mass planets in multiple-planet systems migrated in- 
wards then it seems that we have (at least) two types of pro- 
cesses which may be of importance in planet migration. One 
type of process changes the obliquity while the other does 
not. We could identify disk migration with the latter, and dy- 
namical interactions with the former. We note, though, that 
it is not necessarily true that the compact multiplanet systems 
have experienced inward migration; see, for example, Chiang 
&Laughlin 2012. 

A number of mechanisms have been proposed for changing 
the orbital inclination of a planet. Two processes which have 
attracted particular attention are planet-planet scattering (e.g. 
Rasio & Ford 1996; Chatterjee et al. 2008; Matsumura et al. 
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2010; Nagasawa & Ida 2011) and Kozai Cycles with Tidal 
Friction (KCTF, Eggleton & Kiseleva-Eggleton 2001). Kozai 
cycles can be induced by the influence of a distant stellar com- 
panion (e.g. Wu & Murray 2003; Fabrycky & Tremaine 2007; 
Naoz et al. 2012) or by a distant planet (e.g. Naoz et al. 201 1). 
To confirm Kozai migration via a stellar companion, searches 
for stellar companions to hot-Jupiter hosts will be helpful (e.g. 
Narita et al. 2012). Another way to test theories of hot Jupiters 
involving tidal circularization from a highly eccentric orbit 
is to search for their putative high-eccentricity progenitors. 
Socrates et al. (2012) predicted that there should be a stream 
of gas giant on very eccentric orbits (eccentricity < 0.9) if hot 
Jupiters were transported directly inwards from beyond the 
snow line. Dawson et al. (2012) searched the Kepler database 
for such objects, did not find any, and placed an upper bound 
on such a population. This indicates that if KCTF is an im- 
portant migration path, then likely the starting point for KCTF 
is closer than the snow line. This suggests that gas giants mi- 
grate via a combination of processes. For example initial scat- 
tering or disk migration followed by Kozai cycles and finally 
tidal circularization. 



6. SUMMARY 

In the multiple-planet systems KOI-94 and Kepler-25 we 
measured good alignment between the stellar rotation axes 
and the orbital plane of the transiting planets. For both 
systems we used the radial velocity anomaly (Rossiter- 
McLaughlin effect) during planetary transits to determine the 
degree of alignment. For KOI-94 our result is consistent 
with an independent study by Hirano et al. (2012b). As the 
Rossiter-McLaughlin effect in the Kepler-25 system has only 
a small amplitude, we further measured the distortion of the 
stellar absorption lines directly in another part of the obtained 
stellar spectra. We found consistent results with both meth- 
ods. 

Combining our results with measurements in three other 
multiple planet systems (Kepler-30, Kepler-50, Kepler-65) we 



can now compare the obliquity distributions of multiple planet 
systems to the obliquities measured in Hot- Jupiter systems. 
We find that there is only a 0.0003 % chance that multiple 
planet systems have a isotropic obliquity distribution. This 
is in contrast to the apparent isotropic obliquity distribution 
in Hot- Jupiter systems when taking tidal realignment into ac- 
count (i.e., omitting systems with relatively strong tidal inter- 
actions). 

Our results support the idea that the inward migration of 
close in gas-giants is fundamentally different from the migra- 
tion occurring in compact multiple planet systems. It suggests 
that the planets we see in the multiple planet systems might 
have migrated via disk-planet interactions while Hot-Jupiters 
must have taken a different route. Their path not only brought 
them close to their host stars but also transported them out of 
the orbital plane of the disk out in which they have formed. 
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